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ABSTRACT 

We report on the analysis of the first complete far-infrared spectrum (43- 
197/xm) of the Seyfert 2 galaxy NGC 1068 as observed with the Long Wavelength 
Spectrometer (LWS) onboard the Infrared Space Observatory (ISO). In addition 
to the 7 expected ionic fine structure emission lines, the OH rotational lines at 
79, 119 and 163//m were all detected in emission, which is unique among galax- 
ies with full LWS spectra, where the 1 19 pm line, where detected, is always in 
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absorption. The observed line intensities were modelled together with ISO Short 
Wavelength Spectrometer (SWS) and optical and ultraviolet line intensities from 
the literature, considering two independent emission components: the AGN com- 
ponent and the starburst component in the circumnuclear ring of ~ 3kpc in size. 
Using the UV to mid-IR emission line spectrum to constrain the nuclear ionizing 
continuum, we have confirmed previous results: a canonical power-law ionizing 
spectrum is a poorer fit than one with a deep absorption trough, while the pres- 
ence of a big blue bump is ruled out. Based on the instantaneous starburst age 
of 5 Myr constrained by the Br 7 equivalent width in the starburst ring, and 
starburst synthesis models of the mid- and far-infrared fine-structure line emis- 
sion, a low ionization parameter (U=10~ 35 ) and low densities (n=100 cm -3 ) are 
derived. Combining the AGN and starburst components, we succeed in model- 
ing the overall UV to far-IR atomic spectrum of NGC 1068, reproducing the line 
fluxes to within a factor 2.0 on average with a standard deviation of 1.4. The 
OH 119 pm emission indicates that the line is collisionally excited, and arises in 
a warm and dense region. The OH emission has been modeled using spherically 
symmetric, non-local, non-LTE radiative transfer models. The models indicate 
that the bulk of the emission arises from the nuclear region, although some ex- 
tended contribution from the starburst is not ruled out. The OH abundance in 
the nuclear region is expected to be ~ 10~ 5 , characteristic of X-ray dominated 
regions. 

Subject headings: galaxies: individual (NGC 1068) - galaxies: active - galax- 
ies: nuclei - galaxies: Seyfert - galaxies: emission lines - galaxies: starburst - 
infrared: galaxies. 


1. INTRODUCTION 

NGC 1068 is known as the archetypical Seyfert type 2 galaxy. It is nearby, luminous 
(Lir = 2 x lO a L 0 Bland-Hawthorn et al. 1997), and it has been extensively observed and 
studied in detail from X-rays to radio wavelengths. With a measured redshift of z=0.0038 
(Huchra et al. 1999) (corresponding to a distance of D=15.2 Mpc for H 0 =75 km s -1 Mpc -1 ), 
it provides a scale of only ~ 74 pc/". A central nuclear star cluster has an extent of ~0.6" 


'ISO is an ESA project with instruments funded by ESA Member States (especially the PI countries: 
France, Germany, the Netherlands and the United Kingdom) and with the participation of ISAS and NASA. 
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(Thatte et al. 1997) and a 2.3 kpc stellar bar observed in the near-IR (Scoville et al. 1988; 
Thronson et al. 1989) is surrounded by a circumnuclear starburst ring. Telesco et al. (1984) 
found that the infrared emission in NGC 1068 was due to both the Seyfert nucleus (which 
dominates the 10/xm emission) and to the star forming regions in the bright ~ 3kpc cir- 
cumnuclear ring (which emits most of the luminosity at A > 30 pm). A Bry imaging study 
(Davies, Sugai & Ward 1998) showed a similar morphology and indicated that a short burst 
of star formation occurred throughout the circumnuclear ring of 15-16" in radius within the 
last 4-40 Myr. CO interferometer observations revealed molecular gas very close to the nu- 
cleus (~0.2") suggesting the presence of ~ 10 8 M© within the central 25pc (Schinnerer et 
al. 2000). Recent high resolution H 2 line emission mapping indicates the presence of two 
main nuclear emission knots with a velocity difference of 140 km/s, which, if interpreted as 
quasi-keplerian, would imply a central enclosed mass of 10 8 M© (Alloin et al. 2001). 

In this article, we present the first complete far-infrared spectrum from 43 to 197/rm 
showing both atomic and molecular emission lines (§2). We model the far-IR continuum 
emission using a radiative transfer code and gray body functions (§3) and model the com- 
posite UV- to far-IR atomic emission line spectrum, from our data and the literature, using 
photoionization models of both the active nucleus (§4.1) and the starburst components (§4.2- 
4.3). Moreover, two different non-local, non-LTE radiative transfer codes have been used to 
model the OH lines (§5). Our conclusions are then given in §6. 


2. OBSERVATIONS 

NGC 1068 was observed with the Long Wavelength Spectrometer (LWS) (Clegg et al. 
1996) on board the Infrared Space Observatory (ISO) (Kessler et al. 1996), as part of the 
Guaranteed Time Programme of the LWS instrument team. The full low resolution spectrum 
(43-197//m) of NGC 1068 was collected during orbit 605 (July 13, 1997). Two on-source 
full scans (15,730 seconds of total integration time) and two off-source (6’ N) scans of the 
[CII]158/im line (3,390 seconds of total integration time) were obtained. On- and off-source 
scans had the same integration time per spectral step. Because of the design of the LWS 
spectrometer, simultaneously with the 158^m data, a short spectral scan of equal sensitivity 
to the on-source spectrum was obtained at sparsely spaced wavelengths across the LWS 
range. 

The LWS beam is roughly independent of wavelength and equal to about 80 arcsec. 
The spectra were calibrated using Uranus, resulting in an absolute accuracy better than 
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30% (Swinyard et al. 1996). The data analysis has been done with ISAP 2 , starting from 
the auto-analysis results processed through the LWS Version 7-8 pipeline (July 1998). To 
be confident that newer versions of the pipeline and calibration files did not yield different 
results, we have compared our data with the results obtained using pipeline 10.1 (November 
2001) and we did not find significant differences in the line fluxes or the continuum. 

All the full grating scans taken on the on-source position and the two sets of data on 
the off-source position were separately co-added. No signal was detected in the off-source 
coadds. The emission line fluxes were measured with ISAP, which fits polynomials to the 
local continuum and Gaussian profiles to the lines. In all cases the observed line widths 
were consistent with the instrumental resolution of the grating, which was typically 1500 
km/sec. The integrated line fluxes measured independently from data taken in the two scan 
directions agreed very well, to within 10%. The on source LWS spectrum that resulted from 
stitching the ten LWS channels together using small multiplicative corrections in order to 
match the overlapping regions of each channel with its neighbors is shown in Fig. 1. LWS 
spectra of sources that are very extended within the instrument beam or that peak off center 
are typically affected by channel fringing in the continuum baseline (Swinyard et al. 1998). 
Fortunately, these spurious ripples are hardly noticeable in our LWS spectrum, presumably 
because the far-IR continuum is centrally concentrated towards the center of the LWS 80" 
beam. 

Besides the LWS observations, we also use the SWS observations presented by Lutz 
et al. (2000), to extend the wavelength and ionization-level coverage. Table 1 presents all 
the ISO line flux measurements including those from the SWS with their respective aperture 
sizes. 


3. THERMAL CONTINUUM SPECTRUM 

We have modeled the mid- and far-infrared spectral energy distribution (SED) of NGC 1068 
by means of a non-local, spherically symmetric, radiative transfer code (Gonzalez- Alfonso & 
Cernicharo 1997, 1999). We modeled the nuclear continuum emission closely following the 
approach by Cameron et al. (1993), and using the ISO-SWS fluxes reported by Lutz et al. 


2 The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument 
Teams and Data Centers. Contributing institutes are Centre d’Etudo Spatiale des Rayonnements (France), 
Institute d’Astrophysique Spatiale (France), Infrared Processing and Analysis Center (United States), Max- 
Planck-Insitut fiir Extraterrestrische Phvsisk (Germany), Rutherford Appleton Laboratories United King- 
dom) and the Space Research Organization, Netherlands. 
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(2000). A central source with luminosity 1.5 x 10 u L© illuminates a spherically symmetric 
dust envelope with inner and outer radii of 1 and 200 pc, respectively. The envelope is 
divided into a set of spherical shells where the dust temperature is computed assuming that 
heating and cooling are equal. We assumed a standard silicate/ amorphous carbon mixture 
with optical constants given by Draine (1985) and Preibisch et al. (1993). The density 
profile was assumed to be oc r~ 0 , with (3 and the density at the inner radius regarded as free 
parameters. The predicted flux densities were compared with the continuum observed in 
the SWS spectrum by assuming that the bulk of the continuum in the 2-45 /zm range arises 
from the nuclear region (see below). A reasonable fit to the observed mid-infrared SED is 
shown in Fig. 1, where (3 = 1 and, assuming a gas-to-dust mass ratio of 10 2 , the averaged 
H 2 density at the inner radius is < n/( H 2 ) >= 1.7 x 10 3 cm -3 , the radial H 2 column density 
is A(H 2 ) = 2.7 x 10 22 cm -2 , and the total mass is M = 2.1 x 10 7 M©. The values of A^(H 2 ) 
and M are in agreement with those inferred by Heifer & Blitz (1995). 

The main weakness of our nuclear model is that the flux density at 30-50 gtm is somewhat 
overestimated, a problem that could be more severe if there is a contribution to the continuum 
emission in the SWS aperture by the extended starburst at A > 30 gm . On the other hand, 
the decomposition of the mid-infrared emission performed by Le Floc’h et al. (2001) with 
ISOCAM into the nuclear and extended contributions attributes to the nucleus, between 2 
and 16 gm, significantly higher flux densities than those reported by Lutz et al. (2000) within 
the SWS aperture, suggesting that uncertainties in absolute calibration are affecting one or 
both data sets. Although our continuum model for the nuclear region is thus reasonable 
within the uncertainties, the nuclear continuum emission in the 30-50 gm range could be 
overestimated by a factor as high as 2. The latter is relevant for the analysis of the OH 
emission because, as we argue in section §5, the relative contributions of the nuclear and 
the extended components of the 30-50 gim continuum are important in discerning which 
component is responsible for the OH emission lines. 

The complete far-infrared continuum emission is fitted in Fig. 1 by adding two com- 
ponents at 32 K and 20 K (Spinoglio, Andreani & Malkan 2002). These are gray body 
functions with a steep (/? = 2) dust emissivity law. We identify the 32 K component with 
the starburst ring that is well delineated by the PAH 7.7 gm emission (Le Floc’h et al. 2001). 
Assuming a spherical shell with radius of 1.5 kpc and thickness of 0.3 kpc, the average H 2 
molecular density of the 32 K component is ~ 5 cm -3 . The total mass is 1.6 x 10 9 M© and 
2 x 10 9 M© for the 32 K and 20 K components, respectively. These estimates are in good 
agreement with the ~ 4 x 10 9 M© derived by Planesas et al. (1991) for the molecular ring. 
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4. THE FINE STRUCTURE LINES 

To be able to better constrain the modeling of the line emission of NGC 1068, we have 
combined our far-infrared fine structure line measurements (Table 1) with ultraviolet, optical 
and infrared spectroscopic data from the literature (Kriss et al. 1992; Marconi et al. 1996; 
Thompson 1996; Lutz et al. 2000). The complete emission line spectrum of NGC 1068 from 
the ultraviolet to the fax-IR includes several low-ionization lines that are primarily produced 
outside the narrow line region (NLR) of the active nucleus, as well as intermediate ionization 
lines that originate from both staxburst and AGN emission. For this reason, we find that 
no single model satisfactorily explains all the observed emission lines. We identify two main 
components: 


- an AGN component (the NLR), exciting the high ionization lines and contributing 
little to the low-to-intermediate ionization lines; 

- a starburst component in the circumnuclear ring of the galaxy (e.g. Davies, Sugai & 
Ward 1998) that produces the low ionization and neutral forbidden lines and some of 
the emission in the intermediate ionization lines. This component should also produce 
emission associated with photo-dissociation regions (PDRs) (e.g. Kaufman et al. 1999), 
at the interface with the interstellar medium of the galaxy. 

In this section, we will examine separately the two components that produce the total 
fine structure emission line spectrum of NGC 1068, namely the AGN and the starburst, for 
which we propose two different computations, and we add together these components to 
reproduce the overall observed spectrum from the UV to the far-IR in §4.3. 


4.1. Modeling the AGN 

The first photoionzation model predictions of the mid to far-infrared emission line spec- 
tra of the Narrow Line Regions (NLR) of active galaxies were presented by Spinoglio & 
Malkan (1992), well before the ISO observations could be collected. Alexander et al. (2000) 
used the observed high ionization emission lines to model the obscured ionizing AGN contin- 
uum of NGC 1068 and found that the best-fit spectral energy distribution (SED) has a deep 
trough at 4 Rydbergs, which is consistent with an intrinsic “big blue bump” that is partially 
obscured by ~ 6 x 10 19 cm -2 of neutral hydrogen interior to the NLR. Following their results, 
we have simulated their models, although using a different photoionization code, CLOUDY 
(Version 94.00 Ferland 2000), and then we have varied the shape of the ionizing continuum 
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to include the ionizing continuum derived in Pier et al. (1994). Our goal was to test if the 
Alexander et al. (2000) results were unique and to fit the remaining emission by a starburst 
component, and thereby to derive a composite model of the complete emission line spectrum 
of NGC 1068. 

Specifically, we explore three plausible AGN SEDs. Model A assumes the best fit ionizing 
spectrum derived by Alexander et al. (2000), i.e. with a deep trough at 4 Rydberg (log / = 
-27.4, -29.0, -27.4, -28.2 at 2, 4, 8 and 16 Ryd, respectively). An intrinsic nuclear spectrum 
of NGC 1068 has also been inferred by Pier et al. (1994). Model B assumes the original 
ionizing spectrum derived from Pier et al. (1994). Model C assumes an SED with a Big 
Blue Bump superposed on the Pier et al. (1994) ionizing continuum (log / = -25.8, -25.8, 
-25.8, -27.4 at 2, 4, 8 and 16 Ryd, respectively) as expected for the thermal emission of an 
accretion disk around a central black hole. These three AGN ionizing continua are plotted in 
Fig. 2. For each of models A, B, and C, we have used two component models with the same 
parameters as in Alexander et al. (2000): component 1 has a constant hydrogen density of 
10 4 cm -3 , an ionization parameter U=0.1, a covering factor c=0.45, a filling factor of 6.5 x 
10~ 3 with a radial dependence of the form r -2 , and extends from ~ 21 to ~ 109 pc from 
the center; component 2 has a density of 2 x 10 3 cm -3 , an ionization parameter U=0.01, 
a covering factor of c=0.29, a filling factor of 6.5 x 10~ 4 without any radial dependence, 
and extends from ~ 153 to ~ 362 pc from the center. We have also assumed the “low 
oxygen” abundances adopted by Alexander et al. (2000), in order to be able to compare our 
results with theirs. The inner and outer radii of the emission regions of the two components, 
21, 109, 153 and 362 pc. correspond to angular distances of about 0.26, 1.4, 1.9 and 4.5 
", respectively. Table 2 reports the predicted line fluxes of the three AGN models, A, B, 
and C, together with the observed line fluxes: the line fluxes are given for each of the two 
components 1 and 2, which are treated as independent, and the total flux for each model is 
simply the sum of the fluxes of the two components. 

We can see from Table 2 that only the AGN A and B models, and not the AGN C 
model, reproduce most of t he observed high ionization line fluxes. The low and intermediate 
ionization lines, are expected to have partial or full contributions from starburst and PDR 
components (see §4.2). This first result rules out the presence of a “big blue bump” in 
the ionizing continuum of NGC 1068. To be able to compare the modeled ultraviolet and 
optical lines with the observations, we also listed in Table 2 and 3 their dereddened fluxes, 
assuming two values for the extinction: E b-v — 0.4 mag (Malkan & Oke 1983) and Eg_y 
= 0.2 mag (Marconi et al. 1996). We find that the AGN B model overpredicts several 
of the intermediate ionization lines, such as [SIV]10.5/jm, [NeIII]15.6/rm and [SIII]18.7/um, 
and this discrepancy increases when adding the starburst component because these lines 
are also copiously produced by that component (see next section). On the other hand, the 
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[NeII]12.8/im emission is underpredicted so much so that even with the inclusion of the 
star burst component it cannot be reproduced with this model. As we discuss further in §4.3, 
a composite AGN/starburst model using AGN model A is able reproduce the [NeII]12.8pm 
emission to within the measurement errors. 


4.2. Modeling the starburst ring 

NGC 1068 is known to emit strong starburst emission from the ring-like structure at 
a radial distance of 15 — 16" from the nucleus (total size of ~ 3kpc), traced for example 
by the Br 7 emission (Davies, Sugai & Ward 1998). Mid-IR line imaging observations 
of NGC 1068 have been published by Le Floc’h et al. (2001) based on ISOCAM CVF 
observations. They presented an image of the 7.7 pm PAH feature that shows constant 
surface brightness above the 4th coutour near the nucleus. This suggests that star formation 
is oceuring in the direction of the nucleus so that nuclear spectra will include some emission 
from star formation. In the case of the SWS observations that we are modeling (reported by 
Lutz et al. (2000)), three apertures were used at different wavelengths with the two largest 
also including portions of the brighter starburst ring (see Table 1). To estimate how much 
of the starburst emission is contained in the different apertures used in the observations, 
we have used a continuum subtracted image in the 6.2 pm feature produced by C. Dudley 
(private communication) using the same ISOCAM CVF data set examined by Le Floc’h et 
al. (2001). The 6.2 pm feature is more isolated than the 7.7 pm feature, which is blended 
with the 8.7 pm feature and the silicate absorption feature, but the image compares well with 
the published 7.7 pm image though we have zeroed out the residuals in a 3x12 arcsec 2 region 
centered on the nucleus. Based on this image, the SWS 14x20. 14x27 and 20x33 arcsec 2 
slits contain 13, 23 and 46/i of the 6.2 pm flux contained in the LWS beam respectively, 
without correction for the neglected region of poor residuals (oriented at 45° to our synthetic 
SWS slits). Since PAH features are thought to be a good tracer of PDRs and their associated 
startbursts, we adopt these percentages in our model predictions of SWS line strengths in 
the starburst models presented in this section. Further, in fitting our starburst models to 
the observations, we have abandoned the requirement to obtain the absolute fluxes observed 
at earth, but rather normalized the emission to the intensity of the [Nil] 122pm line observed 
with the LWS, because it has virtually no emission from any AGN model. 

We have chosen the starburst synthesis modeling program Starburst99 (Leitherer et 
al. 1999), to produce input ionizing spectra for the CLOUDY photoionization code. We 
compared the predictions of two different ionizing continua shown in Fig. 3. The first is 
based on an instantaneous star formation law, with a total mass of M = 10 6 M©, while the 
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second is based on a continuous star formation law with a star formation rate of 1 M© yr -1 . 
For both types of models we adopted an age of 5 Myr, a Salpeter IMF (a=2.35), a lower 
cut-off mass of 1 M©, an upper cut-off mass of 100 M©, solar abundances (Z=0.020) and 
nebular emission included. These particular ionizing continua were selected because they are 
consistent with the Br 7 equivalent width observed by Davies, Sugai & Ward (1998) in the 
starburst ring. We have estimated that the Br 7 equivalent width in each of the individual 
regions of the map of Davies, Sugai &; Ward (1998) is in the range 110-180 A. According to 
the Leitherer et al. (1999) models (see their figures 89 and 90), for a value of log(W(Br 7, 
A)) > 2 only instantaneous models with ages less than ~ 6x 10 6 yrs are allowed; similarly, 
only continuous models with ages less than ~ 10 7 yrs are permitted. 

We report in Table 3 the line fluxes predicted for six different starburst models, chosing 
the above instantaneous star formation law model as the input ionizing continuum and using 
CLOUDY with densities of n# = 10, 100, 1000 cm -3 and ionization parameters of Log U 
= -2.5, -3.5. We have also run models with the continuous star formation law presented 
above, but we do not list their results in Table 3, because the differences in the line flux 
predictions are insignificant, while line fluxes strongly depend on the choice of density and 
ionization parameter, as can be seen from Table 3. This result is not surprising because the 
two ionizing continua as derived from the starburst models are in fact quite similar in their 
shape and we are not computing absolute line intensities, but only their ratios to a particular 
observed line. We have also tried continuous starburst models with much longer ages (10, 
20 and 100 x 10 6 years) but, because the shape of the ionizing continuum again does not 
change significantly, the resulting emission line spectrum was indistinguishable from that 
one derived from the models with an age of 5 x 10 6 years. 

In all models the abundances were those of HII regions with grains included. The 
integration was allowed to run until the temperature of the gas in the cloud cooled to T=50 
K in order to include the photodissociation regions present at the interfaces of HII regions 
and molecular clouds. 

It is clear from Table 3 that the models with the higher ionization parameter (log U = 
-2.5) can easily be ruled out, because their emission in many intermediate ionization lines is 
far too high (see e.g. [OIV]26/jm, [OIII]51 ,88/j.m, [NIII]57/im). Moreover, the intermediate 
density of n h = 100 cm -3 gives the best fit to the observed lines, taking into account that 
the AGN component must be added to reproduce the total flux as shown §4.3. 

We estimate the average PDR parameters using the models of Kaufman et al. (1999) 
and the contour plots in Luhman et al. (2003), the measured [C II] 158 and [O I] 145 ^im line 
fluxes (but not the [O I]63 /im line flux which may be affected by absorption and/or shocks), 
and the FIR flux integrated over the LWS spectrum, which we find to be 1.3 x 10~ 8 ergs 
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cm -2 sec -1 . Here we assume that the [C II] line emerges predominantly from PDRs due 
to the strong starburst, rather than the diffuse ionized medium. With this assumption, the 
average PDR gas density and UV radiation field are n h 2 ~ 1000 and Go ~ 300 respectively. 
We note that if instead we assume that the [C II] line flux is dominated by the diffuse ionized 
medium, using the correction factor estimated by Malhotra et al. (2001), we obtain a similar 
gas density nn 2 ~ 1500 but a significantly higher interstellar radiation field Go ~ 1500. For 
both cases, the parameters derived are in the range of those of the normal galaxies in the 
Malhotra et al. (2001) sample, consistent with the assumption that most of the FIR flux 
originates in the starburst ring. 


4.3. Adding the two components 

Summing the line intensities of each one of the two components, the composite spectrum 
of NGC 1068 can be derived and compared with the observed one. We have chosen three 
combinations to compute the composite models, each one with a different AGN model, while 
we adopted the starbust model with n h = 100 cm -3 and Log U = -3.5: 1) the first one (that 
we name CM1, for Composite Model 1) with the AGN ionizing continuum as suggested by 
Alexander et al. (2000) (model AGN A); 2) the second (CM2) with the original Pier et al. 
(1994) (model AGN B); 3) the third (CM3) with the hypothetical bump (model AGN C). 
The results of these three composite models are given in Table 4, compared to the observed 
and dereddened values, assuming the two choices for the extinction (see §4.1). We also show 
the results of the three composite models in a graphical way in Fig. 4, where the modeled to 
the observed flux ratio is given for each line for the case of an extinction of E£_y=0.2 mag. 

A simple x squared test of the three models resulted in a reduced x square of 11.9, 24.6 
and 241 for the three models CM1, CM2 and CM3, respectively. We note that for CM1 the 
line fluxes are reproduced to within a factor of 2.0 on average with a standard deviation of 
1.4 including all of the detected lines for the model with Eb-v=0-2 mag. The infrared fines 
only are reproduced to within a factor of 1.6 on average with a standard deviation of 0.9. 


5. THE OH LINES 
5.1. General remarks 

In NGC 1068, we detect three of the OH rotational lines, all in emission. As shown in 
the energy level diagram of Fig. 5, two of them are fundamental fines, connecting the ground 
state 2 n 3/2 3/2 level with the 2 II 3 / 2 5/2 (the in-ladder 119 pm fine) and with the 2 I1 i/ 2 1/2 
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level (the cross ladder 79 //m line). The third line is the lowest transition of the 2 n 1 / 2 ladder: 
the 163//m line between the J=3/2 and J=l/2 levels. The detected line fluxes are given in 
Table 1. The fact that these three lines are all in emission is in striking contrast with the 
OH lines observed in other bright infrared galaxies, such as Arp 220 (Fischer et al. 1999; 
Gonzalez- Alfonso et al. 2004), Mrk 231 (Harvey et al. 1999), NGC 253 (Bradford et al. 1999), 
and M 82 (Colbert et al. 1999), in which the 119 //m fundamental is in absorption. The 79/xm 
line is sometimes seen in emission and sometimes in absorption; the 163/zm line is always 
seen in emission. In addition to the detections, the ISO-LWS and SWS observations provide 
upper limits on fluxes of the other four lines that arise between the lowest six rotational 
levels. The LWS spectra in the vicinity of the detected lines (and of one of the upper limits), 
are shown in detail in Fig. 6 (histograms). In this section we discuss the physical conditions 
necessary to excite these lines, their probable location within NGC 1068, and detailed model 
fits to their fluxes. A comparison between the observed and modeled line fluxes is given in 
Table 5 and shown in Fig. 6. 


5.2. The excitation mechanism of the OH lines 

The unique OH emission line spectrum of NGC 1068 can provide a powerful way to help 
discriminate between the properties of the molecular clouds in NGC 1068 and the clouds 
in other galaxies in which OH has been observed. Before describing our detailed radiative 
transfer calculations, it is instructive to discuss some conclusions that are model-independent. 
The emission in the OH n :{ / 2 5/2 — 3/2 line at 119 //m cannot be explained by absorption of 
far-infrared photons followed by cascade down to the upper n 3 / 2 5/2 level of the transition. 
Rather, we argue that collisional excitation dominates. Figure 5 shows the energy level 
diagram of OH. There are only two possible paths to excite the 119 gm line via absorption 
of far-infrared photons: via the 35 gm and/or the 53 //in ground-state lines. Excitation by 
either of these routes has other observable consequences. In the case of simple radiative 
cascading, the Einstein-.-! coefficients of the lines involved in the corresponding cascades are 
such that if the 35 /zm absorption path were responsible for the observed 119/zm line flux, 
then the OH Hi/ 2 5/2 — 3/2 line at 98.7 //m would be approximately 5 times stronger than 
the 119 /zm line, while the 98.7 /zm line is not detected. Hence this possibility is ruled out. 
Similarly, if absorption in the 53/zm line were responsible for the observed 119/zm line flux, 
then the 163/zm line would be about 5 times stronger than the 119 /zm line, which it is 
not. We can therefore conclude from the constraints provided by the other far-infrared OH 
lines that the 119/zm emission line is not the result of radiative absorption and cascading. 
The implication is that OH excitation through collisions is more important in NGC 1068 
than in the other observed galaxies and therefore that the gas responsible for the observed 
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emission in the 119pm line resides predominantly in relatively dense and warm enviroments 
in comparison with these other sources. 

The other two observed emission lines, unlike the 119 pm line, need not be collisionally 
dominated. In the case of the II 1/2 3/2— 1/2 163pm line, the most likely excitation mechanism 
is absorption of photons emitted by dust in the 53 and 35 pm lines followed by radiative 
cascade. The upper level of this transition is 270 K above the ground state (Fig. 5), so 
that excitation through collisions is expected to be ineffective in this line. The excitation 
mechanism of the II 1/0 — n 3( / 2 1/2 — 3/2 79 pm line could be a mixure of collisional and 
radiative pumping. The upper level of this transition is 182 K above the ground state, so 
that a warm and dense region could, at least partially, excite the line through collisions. 
Nevertheless, the line could be also excited through the same infrared pumping mechanism 
that results in the observed 163 pm line emission. 

In conclusion, the 119 pm line is collisionally excited, whereas absorption of photons 
emitted by dust in the 53 and 35 pm lines probably dominates the excitation of the 163 pm 
line. The 79 pm OH line may in principle be excited through both mechanisms. 


5.3. Constraints on the spatial origin of the 119 pm OH line 

In NGC 1068 two regions with very different physical conditions can account for the 
observed OH emission as discussed above: the compact nuclear region, and the ring and bar 
where the starburst is taking place. A warm and dense region is required to account for the 
observed 119 pm line emission, given that the line is collisionally excited, so the warm and 
dense nuclear region around the AGN should be considered a good candidate, despite its 
small size (~ 5"; e.g. Planesas et al. 1991; Schinnerer et al. 2000), for the following reasons: 

(i) It is warm: there are ~ 10 3 M 0 of hot H 2 (~ 2000 K) distributed over ~ 5" (Blietz et al. 
1994). From the CO (4-3) to (1-0) line intensity ratio, Tacconi et al. (1994) derive ~ 80 
K for the bulk of the molecular gas, with a mass of ~ 3 x 10 7 M Q enclosed in within 
the central 4" (Heifer Sz Blitz 1995). Lutz et al. (2000) have reported the detection of 
pure H 2 rotational lines within the ISO-SWS aperture, and estimated ~ 2.5 x 10 7 M e 
at ~ 200 K, but these lines may also arise, at least partially, from the inner regions of 
the 3 kpc starburst ring. 

(ii) The molecular clouds within the nuclear region are dense, although there is some 
dispersion in the values derived by several authors based on HCN emission: Tacconi et 
al. (1994) derived an H 2 density of ~ 10 5 cm -3 , whereas subsequent observations and 
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analysis by Heifer & Blitz (1995) yielded a density of ~ 4 x 10 6 cm -3 . An intermediate 
density of ~ 5 x 10 5 cm -3 from HCN and CS, and lower for other tracers, has been 
recently derived by Usero et al. (2004). 

(Hi) The OH abundance is expected to reach high values in regions exposed to strong 
incident UV fields (PDRs, Sternberg & Dalgarno 1995), and in particular in X-ray 
dominated regions (XDRs, Lepp & Dalgarno 1996). The remarkable chemistry found 
by Usero et al. (2004) in the circumnuclear disk of NGC 1068 is indicative of an overall 
XDR and suggests a high OH abundance in the nuclear region. 


Given that the OH 119/zm line is collisionally excited, the possibility that the line might 
arise from the nuclear region can be checked by computing the amount of warm gas required 
to account for the observed emission: 


M w ( M 0 ) = 1.6 x 10 7 x 


10- 5 

Y 

'5 x 10 5 cm -3 ' 

v 

4.3 x 10 11 cm 3 s 1 

_X(OH)_ 

A 

L n(H 2 ) J 

A 

<c iu > 


( 1 ) 


where X(OH) is the OH abundance relative to H 2 , and < q u > is the collisional rate for 
excitation from the ground n 3 /2 3/2 level to the n 3 / 2 5/2 one. Equation 1 assumes that, 
although the line could be optically thick, it is effectively optically thin, and makes use of 
the observed flux of 1.2 x 10 -12 erg s -1 cm -2 . The reference value for the collisional rate, 
< Qu >= 4.3 x 10~ n cm 3 s -1 , corresponds to gas at 80 K (Offer et al. 1994); it decreases 
by a factor of « 2.7 for gas at 50 K and increases by a factor of 3 for gas at 200 K. 


The reference OH abundance we use in this estimate, 10 -5 , is the result of two separate 
studies: first, calculations of molecular abundances by Lepp &; Dalgarno (1996) have shown 
that the OH abundance in XDRs is expected to be about two orders of magnitude higher 
than the abundance of HCN and HCO + . The authors in fact suggested the possibility that 
the high HCN/CO ratio observed in the nuclear region of NGC 1068 could be a consequence 
of enhanced X-ray ionization. Second, the possibility of a chemistry dominated by X-rays 
has found support from observations by Usero et al. (2004), who derive abundance ratios 
of HCN, HCO + , and CN in general agreement with predictions for XDRs. Since the HCN 
abundance derived by Usero et al. (2004) is ~ 10 -7 , AT(OH) in XDRs could attain values 
as high as 10 -5 . On the other hand, the density of 5 x 10 5 cm -3 derived by Usero et al. 
(2004) has been adopted as the reference value in eq. 1. Finally, from the continuum models 
discussed in §3, the mass of gas we derive for the nuclear region is ~ 2 x 10 7 M©, similar 
to the value required in eq. 1. From these estimates we conclude that, if the OH abundance 
is as high as ~ 10 -5 (i.e. if the predictions for XDRs are applicable to the nuclear region 
of NGC 1068), the bulk of the OH 119 /rm line could arise there. This possibility would 
naturally explain why NGC 1068 is unique in its 119 ^m line emission among galaxies with 
full LWS spectra. 
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Finally we ask whether the OH 119/mi line could arise from an even more compact 
region, i.e., from a torus with a spatial scale of 1 pc surrounding the central AGN. According 
to typical parameters given by Krolik & Lepp (1989), a torus is expected to be hot (~ 10 3 
K), could have densities of 10 7 cm -3 , and therefore a mass of ~ 10 5 M©. Also, the OH 
abundance is expected to be very high, 5 x 10 -5 — 10~ 4 . Eq. 1 shows that the relatively 
low mass of the torus (about 2 orders of magnitude lower than the entire nucleus) could be 
compensated by the the higher density, OH abundance, and temperature expected there, so 
that this possibility cannot be neglected. 

The reference values for the nuclear abundance and density adopted in eq. 1 are rather 
uncertain (and possibly extreme). The continuum models of §3 indicate that the mass 
associated with the 32 K dust component, which is identified with the starburst ring, is 
1.6 x 10 9 M©. If ~ 5% of this mass corresponds to warm molecular gas rich in OH, the 
amount of extended warm gas is ~ 8 x 10 7 M©. According to eq. 1, the OH emission at 119 
H m can then also be explained as arising in the ring if the associated PDRs, with assumed 
OH abundance of 2 x 10 -6 (Sternberg & Dalgarno 1995; Goicoechea &: Cernicharo 2002; 
Gonzalez-Alfonso et al. 2004), have densities of a few x 10 5 cm -3 . Since Papadopoulos & 
Seaquist (1999b) found that most of the extended molecular gas resides in dense, compact 
clouds, this scenario seems also possible. However, the continuum from the starburst at 119 
fim is strong, so that one expects that eq. 1 is in this case underestimating M w , and the 
quoted physical parameters, X(OH) and n(H 2 ), are lower limits. The effect of dust emission 
is discussed in detail below. 

In conclusion, a definitive answer to the issue of the spatial origin of the OH 119 gm 
emission cannot be inferred from only the flux observed in the 119 //m line. Nevertheless, 
useful constraints on this subject are given: the line could be either explained as arising from 
the nucleus, with a required OH abundance ~ 10~ 5 , or from the extended ring, with OH 
abundance > 2 x 10~ 6 and density > a few x 10 5 cm -3 . Nevertheless, the radiative transfer 
models described below, which take into account the effect of the continuum emission and 
the excitation of the 79 and 163 fim lines, point towards a nuclear origin of the OH emission. 


5.4. Outline of the models 

Analysis of the OH 79 and 163 /im lines requires the use of detailed radiative transfer 
calculations since, as pointed out above, the emission in these lines is expected to be strongly 
influenced by absorption of far-infrared continuum photons. We therefore proceeded to model 
the OH lines with two different codes, and confirmed that the results were in good agreement 
with each other. One of them, described in Gonzalez-Alfonso & Cernicharo (1997, 1999), 
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has been recently used to model the far-infrared spectrum of Arp 220 (Gonzalez-Alfonso 
et al. 2004), and the other is a Monte Carlo radiative transfer code used as part of a 
detailed study of all the OH lines observed by ISO in galaxies (Smith 2004; Smith et al. 
2004). The latter was developed by the Submillimeter Wave Astronomy Satellite (SWAS) 
mission (Ashby et al. 2000), and is a modification of the original Bernes code. Both methods 
are non-local, non-LTE, assume spherical symmetry, and include a treatment of continuum 
photons from dust mixed in with the gas. Also, both codes take input as a series of concentric 
shells, each of which is assigned a size, gas and dust temperature, H 2 density, velocity and 
turbulent velocity width, and molecular abundance relative to H 2 . The statistical equilibrium 
populations of OH in each spherical shell are computed by including the excitation by dust 
emission, excitation through collisions, and effects of line trapping. We ran two models to 
simulate the nucleus of the galaxy and the starburst extended ring, described in §5.4.1 and 
5.4.2. 


5-4-1 ■ Models for the nuclear emission: constraints on the spatial origin of the 79 and 163 

pm OH lines 

We present models for the nuclear OH emission that implicitly assume that the 119 pm 
emission line arises from the nuclear region: X(OH)= 10~ D is adopted, as well as densities 
~ 5 x 10 5 cm -3 for the bulk of the emitting gas. By assuming a pure nuclear origin for the 
119 pm line, we check whether the other two OH lines could, in such a case, arise from the 
same nuclear region or require a more extended spatial origin. 

The models use dust parameters derived from the nuclear dust model described in §3. 
The gas temperature is assumed to be equal to the dust temperature (Fig. lb). In order to 
avoid overestimation of the continuum flux at the important wavelengths of 35 and 53 pm 
(see section 3), we have decreased the dust abundance by a factor of 2, so that the predicted 
flux in the 35-53 pm range is 40- 50 Jy. 

The H 2 densities derived from the dust model (i.e. a peak density of < n/(H 2 ) >= 
1.7 x 10 3 cm -3 at the inner radius) are not compatible with the densities inferred from 
different molecular tracers. This indicates that the medium is extremely clumped, as has 
been also argued elsewhere (e.g. Cameron et al. 1993; Tacconi et al. 1994). In order to 
account approximately for this dumpiness in our models, the following strategy is adopted: 
we use the “real” n(H 2 ) ~ 5 x 10 5 cm -3 values for the bulk of the gas, and compute the 
volume filling factor /„ =< n(H 2 ) > /n(H 2 ), where the average value is that inferred from 
the dust model. The expected abundances of OH and the dust relative to H 2 , X(OH) and 
AT(dust), are then multiplied by f v , so that the right OH and dust column densities are used 
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in the calculations together with the right density values. The same density profile r 1 that 
was used in the dust model is adopted, so that f v is uniform throughout the nuclear region. 

The modeled fluxes are convolved with the ISO-LWS grating resolution and are com- 
pared with the data in Fig. 6. Solid black lines show the results for the nuclear model that 
assumes X(OH)=10~ 5 and f v — 2.3 x 10~ 5 , the later value implying a density in the outer 
regions (where the bulk of the emission is generated) of 5 x 10 5 cm -3 . The OH column 
density is A’ (OH) = 2.7 x 10 17 cm -2 . Besides the 119 /xm line, the model reproduces the 
emission in the 79 and 163 /xm lines and is consistent with the upper limits given in Table 5. 

We also checked the excitation mechanism of the other OH lines by generating an 
additional model with the same parameters as above except for the continuum emission, 
which is now turned off. In this model, therefore, the lines are excited exclusively through 
collisions with H 2 - The resulting flux of the 119 /xm line remains unchanged when the dust 
emission is ignored, confirming that the line is collisionally excited. On the other hand, the 
flux densities of the 79 and 163 /xm lines decrease in the “pure-collisional” model by factors 
of 2 and 5, respectively, showing that the emission in these lines is much more affected by 
radiative pumping. We conclude that, if the OH abundance in the nucleus were high enough 
to account for the collisionally excited 119 /xm line, the observed fluxes in the 79 and 163 
/xm lines can also be explained as arising in the same nuclear region. 


5-4-2. Models for the starburst emission 

Two simple different approaches have been used to model the OH emission from the 
starburst. First, we have roughly modeled the whole starburst region as a spherical shell 
with external radius of 1.5 kpc, thickness of 0.3 kpc, and average H 2 density < n/(H 2 ) >= 
4.8 cm -3 , so that the corresponding continuum emission is reproduced with T d — 32 K 
(section 3). As shown in section 5.3, the OH 119 /xm emission requires densities of a few x 10 5 
cm -3 , so that we have assumed a volume filling factor f v — 10 -5 and therefore a “real” 
density n(H 2 ) = 4.8 x 10 5 cm" 3 . The kinetic temperature is assumed to be T* — 100 
K, and the OH abundance is X(OH)= 2 x 10“ 6 x f h , where fh — 0.05 is the assumed 
fraction of warm gas. The result of this model is shown in Fig. 6 (upper dotted lines). The 
119 /xm line is reproduced, but the flux densities of both the 79 /xm and 163 /xm lines are 
strongly underestimated. The reason is that the model implicitly assumes that the continuum 
emission, responsible for the excitation of those lines, arises from a very large volume, so 
that the radiation density is weak and has negligible effect on the line excitation. 

Since the OH emission should arise from compact, discrete PDRs in the vicinity of O or 
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early B stars, where the continuum infrared radiation density is expected to be stronger than 
assumed above, we have also tried an alternative approach, which consists of modelling an 
individual “typical” cloud of the starburst. We first model the continuum from an individual 
cloud by assuming a central heating source and computing the dust equilibrium temperatures 
at each radial position that result from the balance of heating and cooling. The continuum 
model is adopted if ( i ) leaving aside a scaling factor (N c , the number of clouds in the 
ensemble), the resulting SED is similar to that of the far-infrared emission of the starburst 
(i.e. the 32 K component found in section 3); the value of N c is determined by requiring that 
the absolute continuum flux from the ensemble of clouds is equal to that observed for the 32 
K component; and (ii) we require that the total mass of the ensemble does not exceed the 
mass inferred from the non-nuclear region (< 4 x 10 9 M©). Once the continuum is fitted, 
calculations for OH are performed by assuming 7* = 7^, and X(OH)= 2 x 10 -6 . 

Several models with various density profiles were found to match the above two require- 
ments. The common characteristic of all of them is the relatively high column density of 
the individual clouds, N(H 2 ) > 10 23 cm -2 , which is a consequence of the low effective dust 
temperature (32 K) of the SED. The results of the simplest model, characterized by a flat 
density profile, n(H 2 ) = 5 x 10 5 cm -3 , are given here for reference. With a radius of 5 x 10 17 
cm, a stellar luminosity of 2 x 10 4 L©, and N c = 4 x 10 6 , the resulting SED is similar to 
that of the 32 K component (§3). For these clouds we obtain a total mass of 1.8 x 10 9 M©. 
The predicted OH emission/absorption is shown in Fig. 6 (lower dashed lines). In spite of 
the relatively high density and temperatures (28-250 K) throughout the cloud, the 119 /xm 
line is predicted to be too weak, and the 79 /xm line is predicted in absorption. We have 
found this result quite general: in models where the OH abundance is high enough and the 
radiation density becomes strong enough to pump the 163 /xm emission, the continuum at 
79 and 119 /xm is absorbed by OH and the predicted emission in the corresponding lines 
is reduced. Models that assume a density profile of r -1 generally predict the 119 /xm line 
in absorption. In some models where the OH abundance was allowed to vary with radial 
position, the 79 /xm line was predicted in emission but by far too weak to account for the 
observed flux density. 

In conclusion, no starburst model is found to reproduce satisfactorily the emission ob- 
served in the three OH lines. If the local infrared radiation density is strong enough to 
pump the 163 /xm line, the other two OH lines are expected to be weak or in absorption. 
Furthermore, the high density assumed for the starburst region would produce a relatively 
high HCN/CO intensity ratio, which is on the contrary ~ 0.01 in the spiral arms (Heifer & 
Blitz 1995). Finally, the PDR models described in 4.2indicateadensityofl-1.5x 10 3 cm -3 , 
i.e. a density much lower than that required to account for the flux density of the OH 119 
/xm line. Therefore, and despite the simplicity of our models, taken together the analysis of 
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the OH lines and the derived PDR parameters indicate that the bulk of the OH emission 
arises from X-ray dominated nuclear regions. 


6. CONCLUSIONS 

The main results of this article can be summarized as follows: 

• The complete far-infrared (50- 200pm) spectrum of NGC 1068 has been observed for 
the first time. The 50-200pm continuum can be interpreted as thermal dust emission 
due to a nuclear component with a luminosity of 1.5 x 10 n L© illuminating a spherically 
symmetric dust envelope with inner and outer radii of 1 and 200 pc, respectively, and 
two additional components at temperatures of 20 and 32K, assuming a steep ( 0 = 2 ) 
dust emissivity law. 

• The far-infrared ISO-LWS spectrum has been complemented with the mid-infrared 
data of ISO-SWS and with shorter wavelength (UV, optical and near-IR) data from 
the literature to assemble a composite atomic spectrum as complete as possible with 
the aim of modeling the different line emission components at work. This approach 
has been necessary especially because of the poor spatial resolution of the ISO spec- 
trometers, which were not able to spatially separate the emission components. The 
lines have been interpreted as arising from two physically distinct components: the 
AGN component and a starburst component, the first one nuclear and the second one 
located in the ring at a radius of 15-16 " from the nucleus. The density and ionization 
parameter of the ~ 5 x 10 h year old starburst are found to be n^ ~ 100 cm -3 and 
log U = -3.5, respectively. Three composite models have been computed with different 
AGN components: the first one has the ionizing continuum as derived from Alexander 
et al. (2000), showing a dee]) trough at energies of a few Rydberg; the second has the 
monotonically decreasing ionizing continuum given by Pier et al. (1994) and the third 
has a “big blue bump” . Two values of the visual extinction ( Eb-v — 0.2 and 0.4) 
have been adopted to correct the optical and ultraviolet line fluxes for the reddening. 
The agreement between the composite model with an AGN ionizing continuum char- 
acterized by the deep trough suggested by Alexander et al. (2000) is very satisfactory, 
taking into account both the simplicity of the phot oionization models chosen to avoid 
dealing with too many free parameters and the large number of lines which originate 
in different physical regimes. The agreement between the observed spectrum and what 
is predicted using the canonical ionizing continuum is poorer, while the presence of a 
big blue bump is ruled out. 
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• The unique OH emission in the 119 /iin line cannot be explained in terms of OH 
excitation through absorption of 35 and 53 /xm photons emitted by dust, but rather 
it is collisionally excited. This indicates the presence of a warm and dense region 
with high OH abundance. A simple excitation analysis yields two main alternatives 
for the spatial origin of the observed 119 /xm line emission: (i) the nuclear region, 
with 2 x 10 7 M© of warm gas (80 K), an average density of n(H 2 ) = 5 x 10 5 cm -3 , 
and an OH abundance of ~ 10 -5 ; (it) the starburst region, if ~ 5% of the associated 
mass (~ 8 x 10 7 M©) is warm (~ 100 K), dense (a few x 10 5 cm -3 ), and rich in OH 
(X(OH)~ 2 x 10 -6 ). 

• Radiative transfer models that simulate the emission/absorption in all the OH lines 
have been performed for both the nuclear and the starburst region. The models for 
the nucleus quantitatively account for the emission in the three OH lines if the nuclear 
physical conditions pointed out above are assumed. On the other hand, no starburst 
model is found to match the three OH lines simultaneously, because the strong far- 
infrared continuum tends to produce absorption, or to weaken the emission, in the OH 
119 and 79 /xm lines (as observed in other galaxies). Therefore, although some con- 
tribution from the extended starburst cannot be ruled out, our models indicate that 
the bulk of the OH emission arises in the nuclear region. The high nuclear OH abun- 
dance required to explain the emission strongly suggest a chemistry deeply influenced 
by X-rays, i.e., an X-ray dominated region. 
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Fig. 1. — a) Spectral energy distribution of NGC 1068 and model fit. Three components, the 
nucleus, the starburst, and a cold component are used, b) and c) Dust temperature versus 
the radial position and radial continuum opacity versus wavelength for the nuclear region. 
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Fig. 2. — The AGN ionizing continua used as input for the photoionization models of 
NGC 1068. The three continua differ in the frequency region between l>E^ yd >100 I while 
outside this region the Pier et al. (1994) spectrum was adopted. The solid line shows the 
continuum derived from Alexander et al. (2000); the dashed line shows a simple power law 
interpolation; the dotted line shows the presence of the predicted ’’big blue bump”. 
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Fig. 3.— The starburst ring ionizing continua used as input for the photoionization models 
of NGC 1068. The two continua are taken from Leitherer et al. (1999) and represent a 
continuous starburst model (solid line) and an instantaneous model (broken line), both with 
5 Myr of age. 
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Fig. 4. — The comparison of the composite models with the observations is shown as the 
ratio of modeled to observed flux ratio for each line, with the ionization potential in the 
x-axis. The assumed reddening is E(B-V)=0.2. Panels from top to bottom: model CM1, 
CM2 and CM3. The short dashed lines represent flux ratios within a factor 3 either ways. 
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Fig. 5. — Energy level diagram of OH. Rotational levels with energies up to 600 K are shown; 
the three lines detected in NGC 1068 are indicated with arrows, as well as the 35 and 53 
fim lines that could play an important role in the radiative excitation. The wavelengths are 
indicated in /zm. A-doubling is ignored because the A-doublets are not resolved with the 
ISO grating resolution. 
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Fig. 6. — Comparison between the observed OH lines and model results. As indicated in 
the upper panel, the upper modeled spectrum (solid lines) corresponds to the model for the 
nucleus, the middle one (dotted lines) corresponds to the starburst modelled as a whole, 
and the lower one (dashed lines) corresponds to the starburst modelled as an ensemble of 
individual clouds (see text for details). 
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Table 1. Measured line fluxes from the LWS and SWS grating spectra, with la 

uncertainties. 


Line 

A 

Flux 

Aperture 

reference 


(pm) 

(10 -13 erg s~ l cm~ 2 ) 

("2) 


[Si IX] 3 P 2 - 3 Pi 

2.584 

3.0 

14 x 20 

1 

[Mg VIII] 2 P 3/2 - 2 P 1/2 

3.028 

11. ± 1.1 

14 x 20 

1 

[Si IX] 3 Pi - 3 Po 

3.936 

5.0 ± 0.6 

14 x 20 

1 

[Mg IV] 2 P 1/2 - 2 P 3/2 

4.487 

7.6 ± 1.5 

14 x 20 

1 

[Ar VI] 2 P 3/2 - 2 P 1/2 

4.529 

15. ± 3. 

14 x 20 

1 

[Fe II] a4 Fg /2 - a6 D 9/2 

5.340 

5.0 

14 x 20 

1 

[Mg VII] 3 P 2 - 3 Pi 

5.503 

13. 

14 x 20 

1 

[Mg V] 3 Pj - 3 P 2 

5.610 

18. ± 2. 

14 x 20 

1 

[Ar II] 2 P 1/2 - 2 P 3/2 

6.985 

13. 

14 x 20 

1 

[Na III] 2 P 1/2 - 2 P 3/2 

7.318 

5.8 

14 x 20 

1 

[Ne VI] 2 P 3/2 - 2 P 1/2 

7.652 

110. ± 11. 

14 x 20 

1 

[Fe VII] 3 P 4 - 3 P 3 

7.815 

3.0 

14 x 20 

1 

[Ar V] 3 P 2 - 3 Pj 

7.902 

< 12. 

14 x 20 

1 

[Na VI] 3 P 2 — 3 Pi 

8.611 

< 16. 

14 x 20 

1 

[Ar III] 3 Pi — 3 P 2 

8.991 

23.0 ± 3.3 

14 x 20 

1 

[Fe VII] 3 P 3 — 3 F 2 

9.527 

4.0 

14 x 20 

1 

[S IV] 2 P 3/2 - 2 P 1/2 

10.510 

58. ± 6. 

14 x 20 

1 

[Ne II] 2 P 3/2 - 2 P 1/2 

12.813 

70. 

14 x 27 

1 

[Ar V] 3 Pi - 3 P 0 

13.102 

< 16. 

14 x 27 

1 

[Ne V] 3 P 2 3 Pi 

14.322 

97. ± 9.7 

14 x 27 

1 

[Ne III] 3 Pi — 3 P 2 

15.555 

160. ± 32. 

14 x 27 

1 

[Fe II] “ 4 P 7/2 - a 4 Fg / 2 

17.936 

< 10. 

14 x 27 

1 

[S III] 3 P 2 — > 3 Pi 

18.713 

40. 

14 x 27 

1 

[Ne V] 3 Pi - 3 P 0 

24.317 

70. ± 7. 

14 x 27 

1 

[O IV] 2 P 3/2 - 2 P V2 

25.890 

190. ± 20. 

14 x 27 

1 

[Fe II] a6 D 7 / 2 - “ 6 D 9/2 

25.988 

8. 

14 x 27 

1 

[S III] 3 Pi - 3 P 0 

33.481 

55. 

20 x 33 

1 

[Si II] 2 P 3/2 - 2 Pi /2 

34.814 

91. 

20 x 33 

1 

[Ne III] 3 P 0 - 3 Pi 

36.013 

18. 

20 x 33 

1 

[O III] 3 P 2 - 3 Pi 

51.814 

114. ± 3. 

80 

2 

[N III] 2 P 3/2 - 2 Pi /2 

57.317 

51.4 ± 2.5 

80 

2 

[O I] 3 Pj — 3 P 2 

63.184 

156. ± 1. 

80 

2 

[O III] 3 Pi 3 P 0 

88.356 

111. ± 1 

80 

2 

[N II] 3 P 2 — 3 Pi 

121.897 

30.5 ± 1.1 

80 

2 

[O I] 3 P 0 - 3 Pi 

145.525 

11.9 ± 0.4 

80 

2 

[C II] 2 P 3/2 - 2 P 1/2 

157.741 

216. ± 1. 

80 

2 

OH 2 ni /2 5/2- 2 n 3/2 3/2 

34.60/34.63 

< 3. 

20 x 33 

2 

oh 2 rii /2 i/2- 2 ni /2 3/2 

79.11/79.18 

14.4 ± 1.5 

80 

2 

OH 2 n 3/2 5/2- 2 n 3/2 3/2 

119.23/119.44 

11.9 ± 1.2 

80 

2 

oh 2 ni /2 3/2- 2 ni /2 i/2 

163.12/163.40 

7.42 ± 0.65 

80 

2 


Note. — (1): from Lutz et al. (2000) and, where errors are available, Alexander et al. (2000); 
(2): this work 
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Table 2. Comparison of observed line fluxes with AGN model predictions 


Line id.A Flux (10“ 13 erg s~ l cm -2 ) 

(fim ) Observed/D 1 /D 2 AGN A model 3 AGN B model 4 AGN C model 5 

Comp. 1 Comp. 2 Comp. 1 Comp. 2 Comp. 1 Comp. 2 


o VI A 1032+.1037 

37.4/4334. /402. 

62.5+33.2 

1. 5+1.4 

12.1+6.4 

0.24+0.20 

155.+80. 

3. 2+1. 7 

(Lya) n A .1215 

101.8/3562./602 

2806. 

1737. 

2951. 

1828. 

3388. 

2182. 

N IV] A .1467 

5.1/103./22.9 

20.0 

4.6 

28.6 

5.0 

67.7 

38.5 

(CIV) n A .1549 

39.7/790. /177. 

488. 

78.5 

462. 

83.2 

1345. 

620. 

Hell A .1640 

21.4/426./95.5 

96.6 

66.5 

96.4 

67.3 

422. 

290. 

[Ne V] A .3426 

15.7/95./38.7 

79.6 

17.0 

56.2 

5.9 

28 8. 

94. 

[Ne III] A .3869+. 3968 

19.2/97./43.2 

38.0+11.4 

29.8+9.0 

145.8+43.9 

88.7+26.7 

162. +49. 

130 +39. 

Hell A .4686 

6.1/27.6/13. 

13.1 

9.6 

13.4 

9.8 

57.8 

41.3 

[O III] A .4959+. 5007 

256 /964-/496 

115.6+334.1 

81.3+235 

290. +837. 

169. +486. 

300. +867. 

229. +661. 

[Si VI] A 1.96 

8.0/9.2/8.6 

10.41 

2.0 

9.7 

2.5 

30.4 

4.8 

[Si VII] A 2.48 

8.3 

9.38 

0.2 

11.8 

.36 

19.0 

0.4 

[Si IX] A 2.584 

3.0 

1.70 


2.6 


1.3 

— 

[Mg VIII] A 3.028 

11. 

6.64 


10.5 


11. 

— 

[Si IX] A 3.936 

5.4 

3.10 


4.9 


2.3 

— 

[Mg IV] A 4.487 

7.6 

3.21 

6.2 

3.5 

5.8 

11.2 

0.5 

[Ar VI] A 4.529 

15. 

14.9 

7.3 

24.1 

6.6 

39.2 

31. 

[Mg VII] A 5.503 

13. 

13.9 

0.23 

15.6 

.14 

33.2 

0.8 

[Mg V] A 5.610 

18. 

5.62 

6.31 

4.0 

3.1 

31.6 

23.9 

[Ar II] A 6.985 

13. 

6.87 

13.4 

1.3 

3.1 

1.4 

3.1 

[Na III) A 7.318 

5.8 

0.66 

0.45 

2.2 

1.4 

1.5 

1.1 

[Ne VI] A 7.652 

110. 

134. 

8.8 

93.7 

2.6 

455. 

54. 

[Fe VII] A 7.815 

3.0 

1.92 

1.90 

1.9 

1.0 

8.6 

6.4 

[Ar V] A 7.902 

< 12. 

2.90 

1.91 

3.3 

2.5 

4.5 

7.0 

[Na VI] A 8.611 

< 16. 

1.10 

0.18 

0.8 

.08 

3.8 

0.8 

[Ar III] + [Mg VII] A 8.991 

25. 

17.8+17. 

9.0+0.36 

18.2+19.4 

19. +.19 

16. +40. 

15.3+1.0 

[Fe VII] A 9.527 

4.0 

7.96 

8.2 

7.7 

4.5 

35.3 

26.8 

[S IV] A 10.510 

58. 

76.5 

43.6 

155. 

117. 

133. 

161. 

[Ne II] A 12.813 

70. 

24.8 

23.9 

2.9 

6.0 

1.9 

4.3 

[Ar V] A 13.102 

< 16. 

2.94 

2.86 

3.4 

3.7 

4.7 

10.3 

[Ne V] A 14.322 

97. 

79.2 

50.7 

74.6 

23.7 

332. 

232. 

[Ne III] A 15.555 

160. 

68.1 

37.6 

136.7 

104. 

91.4 

79.6 

[S III] A 18.713 

40. 

33.7 

41.5 

31.5 

61.9 

28.2 

51.3 

[Ne V] A 24.317 

70. 

34.8 

39.8 

31.2 

18.2 

139. 

182. 

[O IV] A 25.890 

190. 

22.2 

55.4 

25.6 

57.5 

69.6 

200. 

[S III] A 33.481 

55. 

8.59 

26.3 

7.7 

38.1 

6.7 

33.1 

[Si II] A 34.814 

91. 

11.2 

39.5 

8.6 

38.5 

5.9 

37.1 

[Ne III] A 36.013 

18. 

4.97 

3.2 

10.2 

8.8 

6.9 

6.8 

[O III] A 51.814 

no. 

19.2 

27.3 

28.2 

59.5 

19.8 

44.7 

[N III] A 57.317 

51. 

4.64 

8.7 

4.9 

17.4 

3.7 

14.4 

[O I] A 63.184 

156. 

2.50 

3.6 

1.5 

3.7 

0.9 

3.8 

[O III] A 88.356 

110. 

2.53 

6.5 

3.8 

14.2 

2.6 

11.3 

[N II] A 121.897 

30. 

0.13 

0.86 

.06 

.40 

.08 

0.5 

[O I] A 145.525 

12. 

0.16 

0.28 

.09 

.29 

.05 

0.3 

[C II] A 157.741 

220. 

0.46 

3.1 

.29 

2.2 

.11 

1.8 
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Table 2 — Continued 


Line id.A Flux (10 -13 erg s -1 cm" 2 ) 

(/tm ) Observed/D 1 /D 2 AGN A model 3 AGN B model 4 AGN C model 5 

Comp. 1 Comp. 2 Comp. 1 Comp. 2 Comp. 1 Comp. 2 


1 Dereddened line flux, assuming Es_\/ = 0.4 
2 Dereddened line flux, assuming Ej = 0.2 

3 AGN A parameters: component 1: Log U=-l., Log n=4, internal radius ~ 21 pc, external radius ~ 
109 pc, ionizing spectrum from Alexander et al. (2000): component 2: Log U=-2., Log n=3.3, internal 
radius ~ 153 pc, external radius ~ 362 pc, ionizing spectrum from Alexander et al. (2000). 

4 AGN B parameters: same as AGN A models, but with the ionizing spectrum from Pier et al. (1994) 

5 AGN C parameters: same as AGN A models, but with the ionizing spectrum that includes a big blue 
bump (see text) 
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Table 3. Comparison of observed line fluxes with the Ring Starburst model predictions 


Line id. A 
(jan ) 

Observed 

SBR A 1 

Flux (10 
SBR B 2 

13 ergs" 1 
SBR C 3 

cm 2 ) 

SBKD 4 

SBR E 5 

SBR F 6 

O VI A .1032+ 1037 

37.4/4334. /402. 

— 

— 

— 

— 

— 

— 

(Lya) n A .1215 

101.8/3562./602 

221. 

36.3 

284. 

400. 

1135. 

1424. 

N IV] A .1487 

5.1/103./22.9 

22.8 

0.02 

41.8 

0.02 

177. 

0.08 

(CIV) n A .1549 

39.7/790. /177. 

162. 

0.39 

277. 

.75 

1154. 

1.80 

Hell A .1640 

21.4/426./95.5 

235. 

38.9 

373. 

47.8 

1466. 

161. 

[Ne V] A .3426 

15.7/95./38.7 

16.2 

— 

30.9 

— 

125. 

— 

[Ne III] A .3869+. 3968 

19.2/97./43.2 

106.+32. 

24.6+7.4 

165+49.7 

30.4+9.2 

659+199 

105+31.6 

Hell A .4686 

6.1/27.6/13. 

33.5 

5.8 

53.1 

7.1 

207. 

23.8 

[O III] A .4959+. 5007 

256 /964. /496 

439+1269 

19.8+57. 

720+2077 

246+71.2 

2882+8307 

85.8+247. 

[Si VI] A 1.96 

8.0/9.2/8.6 

— 

— 

— 

— 

— 

— 

(Si VII] A 2.48 

8.3 

— 

— 

— 

__ 

— 

— 

[Si IX] A 2.584 

3.0 

— 

— 

— 

— 

— 

— 

[Mg VIII] A 3.028 

n. 

__ 

— 

— 

— 

— 

— 

[Si IX] A 3.936 

5.4 

— 

— 

— 

— 

— 

— 

[Mg IV] A 4.487 

7.6 

1.0 

— 

1.7 

— 

6.6 

— 

[Ar VI] A 4.529 

15. 

1.7 

— 

3.3 

— 

13.3 

— 

[Mg VII] A 5.503 

13. 

— 

— 

— 

— 

— 

— 

[Mg V] A 5.610 

18. 

0.84 

— 

1.6 

— 

6.2 

— 

[Ar II] A 6.985 

13. 

2.6 

3.3 

3.1 

4.0 

11.7 

13.6 

[Na III] A 7.318 

5.8 

0.2 

0.09 

0.3 

0.1 

1.2 

0.4 

[Ne VI] A 7.652 

no. 

4.5 

— 

8.0 

— 

35. 

— 

[Fe VII] A 7.815 

3.0 

0.13 

— 

0.3 

— 

1.2 

— 

[Ar V] A 7.902 

< 12. 

1.9 

— 

4.0 

— 

15.6 

— 

[Na VI] A 8.611 

< 16. 

— 

— 

— 

— 

— 

— 

[Ar III] + [Mg VII] A 8.991 

25. 

16. 

5.6 

23. 

6.9 

89.7 

23.1 

[Fe VII] A 9.527 

4.0 

0.6 

— 

1.1 

— 

4.3 

— 

[S IV] A 10.510 

58. 

91. 

0.6 

157. 

0.8 

589. 

2.4 

[Ne II] A 12.813 

70. 

3.3 

4.9 

3.9 

5.8 

15.2 

20.4 

[Ar V] A 13.102 

< 16. 

5.7 

— 

10.6 

— 

41.4 

— 

[Ne V] A 14.322 

97. 

74.4 

— 

134. 

— 

524. 

— 

[Ne III] A 15.555 

160. 

97.3 

40.3 

143. 

48.7 

552. 

163. 

[S III] A 18.713 

40. 

65.3 

1.6 

99. 

20.3 

386. 

69. 

[Ne V] A 24.317 

70. 

83.6 

— 

148. 

— 

490. 

— 

[O IV] A 25.890 

190. 

496. 

1.4 

828. 

1.7 

2615 

4.7 

[S III] A 33.481 

55. 

199. 

51. 

279. 

58.4 

649. 

125. 

[Si II] A 34.814 

91. 

45.8 

45. 

53. 

49.3 

105. 

100. 

[Ne III] A 36.013 

18. 

17.4 

6.9 

26. 

4.3 

97. 

28. 

[O III] A 51.814 

no. 

921. 

62. 

1559. 

80.4 

4800. 

222. 

[N III] A 57.317 

51. 

349. 

30. 

468. 

31.8 

750. 

44.4 

[O I] A 63.184 

156. 

65. 

99. 

77. 

120. 

300. 

383. 

[O III] A 88.356 

no. 

1539. 

105. 

1812. 

99.5 

1830. 

90. 

[N II] A 121.897 

30. 

30.* 

30* 

30.* 

30* 

30* 

30.* 

[O I] A 145.525 

12. 

6.5 

10. 

7.7 

12. 

27. 

36. 

[C II] A 157.741 

220. 

294. 

377. 

130. 

158. 

1440. 

145. 
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Note. — *: this line was used for normalization 

*SBR A parameters: Log U=-2.5, Log n=1.0, ionizing spectrum from Starburst99 with instantaneous star-formation law, M 
= 10 6 M©, IMF: =2.35 M up = 100M©, M/ ou , = 1 Af©, nebular emission included, Z=0.020, age of 5 Myr. The integration was 
stopped at a temperature of 50K, the adopted abundances are those relative to HII regions and grain emission is included. 

2 SBR B parameters: Log U=-3.5, Log n=1.0, all other parameters as for SBR A. 

3 SBR C parameters: Log U=-2.5, Log n=2.0, all other parameters as for SBR A. 

4 SBR D parameters: Log U=-3.5, Log n=2.0, all other parameters as for SBR A. 

5 SBR E parameters: Log U=-2.5, Log n=3.0, all other parameters as for SBR A. 

6 SBR F parameters: Log U=-3.5, Log n=3.0, all other parameters as for SBR A. 


Appendix C (paper 2) 


-34- 


Table 4. Comparison of observed line fluxes with composite model predictions 


Line id.A 
(/"“ ) 

Flux (10- 13 
Observed/D 1 /D 2 

erg s 1 
CM1 3 

cm 2 ) 
CM2 4 

CM3 5 

O VI A . 1032+1037 

37.4/4334. /402. 

98.6 

18.9 

240. 

(Lya) n A .1215 

101.8/3562./602 

4943. 

5179 

5970. 

N IV] A .1487 

5.1/103./22.9 

24.6 

33.8 

106.2 

(CTV) n A .1549 

39.7/790./177. 

567. 

546 

1967. 

Hell A .1640 

21. 4/426./95.5 

211. 

211.5 

761. 

(Ne V] A .3426 

15.7/95./38.7 

96.6 

62.1 

382. 

[Ne III] A .3869+. 3968 

19.2/97./43.2 

127.8 

344.7 

420. 

Hell A .4686 

6.1/27.6/13. 

29.8 

30.3 

106.2 

[O III] A .4959+. 5007 

256./964./496 

861.8 

1879 

2153. 

[Si VI] A 1.96 

8.0/9.2/8.6 

12.4 

12.2 

35.2 

[Si VII] A 2.48 

8.3 

9.6 

12.2 

19.4 

[Si IX] A 2.584 

3.0 

1.7 

2.6 

1.3 

[Mg VHI] A 3.028 

11. 

6.6 

10.5 

11. 

[Si IX] A 3.936 

5.4 

3.1 

4.9 

2.3 

[Mg IV] A 4.487 

7.6 

9.4 

9.3 

11.7 

[Ar VI] A 4.529 

15. 

22. 

30.7 

70.2 

[Mg VII] A 5.503 

13. 

14. 

15.7 

34. 

[Mg V] A 5.610 

18. 

12. 

7.1 

55.5 

[Ar II] A 6.985 

13. 

24. 

8.4 

8.5 

[Na III] A 7.318 

5.8 

1.1 

3.7 

2.7 

[Ne VI] A 7.652 

110. 

143 

96.3 

509. 

[Fe VII] A 7.815 

3.0 

3.8 

2.9 

15. 

[Ar V] A 7.902 

< 12. 

4.8 

5.8 

11.5 

[Na VI] A 8.611 

< 16. 

1.3 

0.9 

4.6 

[Ar III] + [Mg VII] A 8.991 

25. 

51. 

64.1 

79.2 

[Fe VII] A 9.527 

4.0 

16. 

12.2 

62.1 

[S IV] A 10.510 

58. 

121 

273. 

295. 

[Ne 11] A 12.813 

70. 

55. 

14.7 

12. 

[Ar V] A 13.102 

< 16. 

5.8 

7.1 

15. 

[Ne V] A 14.322 

97. 

130. 

98.3 

564. 

[Ne III] A 15.555 

160. 

155. 

289.4 

220. 

[S III] A 18.713 

40. 

95. 

113.7 

99.8 

[Ne V] A 24.317 

70. 

75. 

49.4 

321. 

[O IV] A 25.890 

190. 

80. 

84.8 

272. 

[S III] A 33.481 

55. 

93. 

104.2 

98.2 

[Si II] A 34.814 

91. 

100. 

96.4 

92.3 

[Ne III] A 36.013 

18. 

12.5 

23.3 

18. 

[O III] A 51.814 

110. 

127 

168.1 

145. 

[N III] A 57.317 

51. 

45. 

54.1 

50. 

[O I] A 63.184 

156. 

126. 

125.2 

125. 

[O III] A 88.356 

110. 

109. 

117.5 

113. 

[N II] A 121.897 

30. 

31. 

30.5 

30.6 

[O I] A 145.525 

12. 

12.4 

12.4 

12.3 

[C II] A 157.741 

220. 

161 

160.4 

160. 


^ereddened line flux, assuming Eb-V = 0.4 
2 Dereddened line flux, assuming Eb~v — 0-2 


Appendix C (paper 2) 





-35- 


3 CM1 = AGN A + SBR D 
4 CM2 = AGN B + SBR D 
5 CM3 = AGN C + SBR D 
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Table 5. Comparison of observed line fluxes with model predictions for the nuclear region 


Line id.A 

Flux (10- 12 
Observed 

erg s 1 cm 2 ) 
Modeled 

Notes 

34 fim 

<0.1 

-0.70 

(absorption) 

53fim 

< 1.2 

-0.84 

(absorption) 

79 fim 

1.44 

1.46 


Mfim 

< 1.2 

0.16 


98fim 

< 1.2 

0.31 


119/xm 

1.19 

1.52 


163 pm 

0.74 

0.51 
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